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The motion on the sphere 5^ in potential V = (xiX2a:3)~^/^ is considered. The Lax represen- 
tation and the hnearisation procedure for this two-dimensional integrable system are discussed. 



1 Description of the model 



The system under consideration is a special case of the following mechanical system in the 
nine- dimensional space 

F*{t) + I det G = . (1.1) 

Here F{t) and G are 3x3 matrices, F* means transpose matrix and 7 is an polytropic index. 
The components Fjj^ of the matrix F are coordinates on the configuration space MP. These 



equations of motion have been studied many authors, see [|1^, |, g, g]. 

According to we shall consider symmetric constant matrices G = G* only. In this 
case, by using canonical transformation of variables F' = UFU*, we can reduce the constant 
matrix G to the diagonal matrix with the following diagonal elements ±1 or 0. Moreover, 
with physical point of view we can put G = I without loss of generality |jl2 



The Newton equations (|1 . 1|) arise in solution of the hydrodynamical equations represent- 
ing dynamics of a cloud of compressible gas expanding freely in an otherwise empty space. 
This model has a rich history related with the names of Dirichlet, Dedekind and Riemann. 
For an extensive discussion of the model we refer to a book which should be taken as a 
general reference on the subject. 

At G = G* the invariance of the problem under rotation and under internal motion of 
the gas leads to conservation of the angular momentum and operator of vorticity 

J = F{t) F*{t) - F{t) F*{t) , K = F*{t) F{t) - F*{t) F{t) . 

Each of J, K is an antisymmetric matrix with 3 independent components. Thus, equations 
( |1.1| ) possess an enlarged symmetry group 5*0(4) ~ 5*0(3) x 50(3). There is also a discrete 
symmetry, which allows the vorticity and the angular momentum to be interchanged. This 
discrete symmetry is identical with the duality principle of Dedekind M . 



For the perfect monatomic gas, at 7 = 5/3 the system of equations (|LT|) possesses one 
more integral of motion 

r2=tr(F*(t)F(t)), at 7 = ^. (1.2) 

The qualitative behavior of solutions at the different values of the adiabatic index 7 and 
some partial cases of motion are discussed in [Q. 

It is easy to separate the diagonal and non-diagonal components of equations ( |1.1|) [0, ||, 
P]. The non-diagonal components give six kinematical equations involving only the inertial 
properties of the gas-cloud but not the pressure force. The diagonal components give three 
dynamical equations determining the rate of expansion of the gas-cloud under the influence 
of the pressure force. 

Below at 7 = 5/3 we shall consider a free expansion of an ellipsoidal gas cloud with the 
fixed orientation, having zero angular momentum and zero vorticity K = J = 0. In this case 
matrix F{t) is diagonal for all t 

F(t) = diag{ Xi, X2, X3 ) (t) 

and the constant matrix G = / is equal to unit. The corresponding three equations of motion 
are given by 

const 5 , , 

xixi = X2X2 = X3X3 = r^TTT, at 7=o- (1-3) 
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Additional integral of motion ( [L.2| ) is equal to the radius of sphere = Yl^l- using 
this integral our system on may be reduced to the system on the sphere S*^. At 7 = 5/3 
the third independent integral of motion was founded in ||^. 

The reduced system has the configuration space diffeomorphic to the Euclidean motion 
group -£^(3) = 5*0(3) x [Q. It allows to identify the phase space of this system on T*S'^ 
with the cotangent bundle T*E{3). The kinetic energy is a left-invariant Riemannian metric 
on E{3). It is determined by some quadratic form on the dual space e*(3) of the Lie algebra 
e(3) 1,0. 

By using the Killing form the dual space e*(3) may be identified with algebra e(3) = 
so{3) © M.^, the semi-direct sum of so(3) and the abelian space M^. Let two vectors J G 
so(3) ^ M.^ and x G be coordinates in the dual space e*(3) equipped with a natural 
Lie-Poisson brackets 

■{t/j , i/j ■} £ijk Jk 1 {Ji 1 ■^j } ^ijk -^k ) 

{xi,Xj} = 0, /c = 1,2,3, . (1.4) 

Here Eijk is the standard totally skew-symmetric tensor. The generic coadjoint orbits of -£^(3) 
in e*(3) are four dimensional symplectic leaves specified by the two Casimir elements 

Ci = {x,x) = XiXi] C2 = {J, x) = JiXi . (1.5) 

Here {x,y) means inner product in R^. Thus, the dual space e*(3) decomposes into the 
coadjoint orbits 

Oc,,c, = {{J,x} G : Ci = ci, C2 = C2} (1.6) 

which are invariant with respect to the usual Euler-Poisson equations in e*(3) |13|. 

Let us introduce a complex analog of the Lie algebra e(3), as a semi-direct sum of so(3, C) 
and the complex space C^. This algebra e(3, C) = so(3, C) © is equipped with the same 
Lie-Poisson brackets ( |1.4|) and the Casimir operators ([L5| ). In contrast with the usual e(3) 



algebra, it allows us to consider non-trivial representations at zero value ci = of the first 
Casimir operator Ci. 

The condition ci = has no obvious physical or geometric meaning. Of course, we can 
not consider real sphere of zero radius, but with mathematical point of view ci is an arbitrary 
value of the Casimir element. Note, that "non-physical" representations of the algebra s/(2) 
with the zero spin s = are useful in physics [0, as well. 



Proposition 1 At zero value Ci = of the first Casimir operator Ci ( \l.di ) the following 
transformation in so* (3, C) C e*(3, C) 

J ^ J = J + 777^^5 OGC (1.7) 

[Xi X2 Xs)^/'^ 

is an outer automorphism of the representation of e(3, C) . 

By using embedding e(3) C e(3, C) let us consider known integrable tops on this com- 
plex algebra e*(3,C). Applying transformation (|1 . 7|) one can get integrable deformations 
of these top on the one-parameter set of orbits Oi (ci = 0,C2 = const). Sometimes outer 
automorphism ( |1.7|) allows us to get much more. 
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As an example, let us consider spherical top with the standard Hamilton function 

H={jJ) = Jl + Jl + Jl (1.8) 

and with the non-standard second integral of motion 

K = JiJ2h (1.9) 

defined on the subalgebra so(3) only. Of course, substituting vector J (|1.7|) one gets inte- 
grable deformation of this symmetric top at Ci = 0. However, we can prove the following 

Proposition 2 Outer automorphism ?| j maps Hamiltonian (Lt_) of the spherical top into 
the following Hamiltonian 

3 

H = Y Jl + 2ia- -a2- (1.10) 

The proposed deformation of the spherical top is completely integrable on the both one- 
parameter sets of orbits 

Oi : (ci = 0, C2 = const) and O2 : (ci = const, C2 = 0) . 

Moreover, to consider the second orbits we can return to the usual real Lie algebra e*(3). 

At Ci = the second integral is the image of known integral K ( |1.9|) 



\Xi X2 X3J 

(1.11) 

+ 7 ^ " M /o (^1 J2 X3 + Ji X2 J3 + Xi J2 J3) -ia^ ■ 

[Xi X2X3)^''^ 

At C2 = the Hamilton function ( |1.10| ) is in involution with the following second integral of 
motion 

K = J,J2J3 + a' (- + - + -) (XlX2X3)^/^ (1.12) 
\Xi X2 X3J 

In contrast with (|1.11|) , here we removed imaginary terms and changed the sign before the 
rest potential term. We do not know origin of such additional transformation as yet. 
In the natural variables 

X 

[Xi X2 X3)^'-^ 

transformation (|1.7|) becomes a shift J = J + iay and the integrals of motion ( |1.10|) and 
( |1.12| ) are given by 

H = Jl + Jl + Jl-a\yl + yl + yl), 

(1.13) 

K = J,J2J3 + a' ( ^ + ^ + ^' 

V yi 2/2 2/3 , 



The Euler-Poisson equations on e*(3) generated by ( p,.10| ) are given by 



d 2 

-^y =-^iy,y) Jxy'\ {y^y) = yl + yl + yl 



:i.i4) 



where x x y means standard vector product in R^. Thus, we rewrite initial very symmet- 
ric equations of motion (|1.3|) defined on the configuration space as the Euler-Poisson 
equations ( |1.14| ) defined on the phase space e*(3). 



2 Lax representation 

The main purpose of this note is to rewrite equations of motion ( |1.14| ) in the Lax form 



|l=[L,M]. (2.15) 

Let us briefiy recall construction of the Lax pair for the Neumann system. The Neumann 
system is an integrable system on the sphere with quadratic potential (see (|L13|) ). Its phase 
space may be modelled on the dual space e*(3) at C2 = 0. The corresponding Euler-Poisson 
equations are equal to 

J = X X z , ■^x = —Jxx, z = —diag{ai, a2, as) X . (2-16) 
at at 

Here aj be arbitrary parameters. 

The Neumann system possesses the necessary number of the quadratic integrals of 
motion. Nevertheless, the Lax pair can not be constructed in framework of the algebra 
e(3) = so(3) © M.^. Namely, for the Neumann system and some others system, we have to 
use the Cartan-type decomposition of the Lie algebra (7/(3, M) = so(3) + Symm{3) |jl3 . 



Let us introduce antisymmetric matrix of angular momentum J7 G so(3) and symmetric 
matrix of coordinates X G Symm{?)) 

J E so{3) - : Jij = EijkJk , 

X G Symm{?)) : Xij = Xi Xj . 

Then the Lax representation for the Neumann system are given by 

L = diag{ai, a2, as) \ + J + X , M = -\-^ X . (2.17) 
Let us present these Lax matrices explicitly 

ai\ \ / -^2 \ 1/^1 X1X2 XiXs 






-J3 


J2 


■h 





-Jl 


-J2 


Jl 






a2X I + I J3 -Jl I + ^ I X1X2 X2 X2X3 
a-sX J \ -J2 Jl / \ xixs X2X3 xl 
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M 



Xi X2 X3 

A 



Xl 



X2X3 
-1 

^3 



X 



-1 



X3 

X2 
Xl X3 

x7' 



X 



^1 



1 



Xl X2 



\ I 

Now let us turn to the deformation of the completely symmetric top ( |1.1CI| ). The Lie 
algebras R'^ with vector product and so(3) with usual commutator may be identified by using 
the Lie algebras isomorphism 



X 




X 






X3 


-X2 


-X3 





Xl 


X2 


—Xl 






G so(3) 



In ( |1.7|) the element J G so(3) has been added with the vector x G M^. Thus, defining outer 
automorphism ( |1.7| ), we implicitly used this property of the three dimensional Euclidean 
space. Below we shall use the same property to construct the Lax representation. 

The main recipe is to rearrange items in decomposition gl{3,M.) = so{3) + Symm{'i). 
Let us introduce symmetric matrix of angular momentum J G Symm{'i) and antisymmetric 
matrix of coordinates X G so(3) 



J G Symm{3) : 



Ji 



1-3 



-'ijk 



Jk : 



X G so(3) 



83 . 



Xi 



^ijk Vk ■ 



where | Sijk \ means absolute value of Sijk- 



Proposition 3 At C2 = the equations of motion l \l-14\ ) on the sphere S*^ generated by the 
Hamilton function dl.lS ) can be written in the Lax form i\2.1L ) with the following matrices 



L = \I + J + aX 

More explicitly, the first matrix is 



Mi. 



2a 



£ijk X, 



















A 





H 




^ 





A j 





J3 J2 
J3 Ji 
J2 Ji 



+ a 






2/3 


-2/2 


-2/3 





2/1 


2/2 


-2/1 






and the second matrix is given by 



M 



2aci 



1/2 



3^(2/, 2/) 



ys' 2/2"' 
' y^' 




2/3_ 

2/2"^ 2/r^ 



(2.18) 



The spectral invariants of L{X) give rise to both integrals of motion in involution 

detL(A) = - HX + 2K . 
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The proposed Lax matrix L{X) has a trivial dependence on spectral parameter A, which is 
similar to the Lax matrix for the Kowalewski top due by Perelomov Therefore, we can 
not construct a suitable spectral curve and can not directly integrate equations of motion. 
Recall, in |]T3|| the Perelomov matrices have been embedded into the general Lax matrices 
with the spectral parameter. It forces us to consider the Lax representation ( p.l8| ) as a first 
attempt to build an adequate Lax pair. We believe the desire Lax pair explains the peculiar 
geometry and the origin of integrability of the considered motion on the sphere S'^. 



3 Linearisation procedure 



In closing this note we briefly discuss results obtained in ||T0[ within the modern theory 
of linearisation of the two-dimensional integrable systems fl], |l6l. The proposed in [1T0| 
procedure of integration has an unusual form, which is closely related with the concrete 
system of equations. This construction is technically closed to the Chaplygin approach to 
the Kirchhoff equations at C2 = 0. On the other hand, modern theory of linearisation allows 
us to consider different integrable systems such as the Neumann problem, Henon-Hailes 
system, Toda lattice, Kowalewski top, Goryachev-Chaplygin top and many other ^, |16 



However, in this common powerful method it is necessary to rewrite equations of motion 
at some suitable variables. These variables have to satisfy the special conditions ^, |TB[ . As 



an example, to integrate the Toda lattice we have to introduce so-called Flaschka variables. 
For other system such variables may be introduced by using the Kowalewski-Painleve analysis 
or the algebro-geometric tools. Nevertheless, if we might have introduced such variables, the 
Adler and van Moerbeke methods ^, |l^ enable us to integrate a given mechanical system. 

The aim of this Section is to introduce an analog of the Flaschka variables for the 
deformations of the spherical top. At these variables we may directly apply the Adler and 
van Moerbeke methods to a given integrable system. These results will be published at the 
consequent publications. 

The three body Toda lattice is the Hamiltonian system defined as 



1 ^ 



Here {pj, qj}^=i be pairs of the canonical physical variables. According to |[T|, |5, 0, in the 
Flaschka variables 

z-^ = e«^-'?2 , Z2 = e'^^-'^' , Z3 = e^^*-"! , Zi ^2 -23 = 1 

Z4 = -pi , Z5 = -P2 , Zq = -P3 

the corresponding equation of motion have the following form 

— = (^5 — Z4) , 

Z2 = Z2 (Zq - Z5) , ^Z5 = Z2-Zi, (3.19) 

at 
d 



^ ^3 — ^3 1^4 — Zq) 



d 






Jt'' 


= Zl 


- Z3 


d 






Jt'' 


= Z2 


- Zl 


d 






dt'' 


= Z3 


- Z2 



7 



The Toda flow has the following four constants of motion 



Zi — Zi Z2 Z3 — 1 , 

Z2 = z^ + Z5 + zq = di = , 

Z-i = \ {zj + zl + zl) + zi + Z2 + zz = ai, (3.20) 

Z4 = ^4 Z^ Zq — Zl Zq — Z2 Z4 — z^ z^ = bi . 

At c?i 7^ the variable Q = qi + q2 + can not be restored from variables {z}^^i. Really, 
we have to add to variables {zj} some other variables with the trivial dynamics [||]. Below 
we shall introduce analog of the Flaschka variables for the integrable deformations of the 
spherical top. 

There is a discrete permutation group acting on the vectors q and p simultaneously 



Vq, pi — >Vp, V = \ 1 \ , V^ = l. (3.21 




According this point symmetry, the invariant manifold defined by (|3.20| ) has the third order 
automorphism given by 

(Zl, Z2, Zs, Z4, ^5, Zq) I > {Z2, Z3, Zl, Zrj, Zq, Z4) . (3.22) 

This automorphism simplifies the Adler and van Moerbeke analysis ^, which applied to 
this system, gives linearisation to the Toda flow. 

Recall, the Kowalewski-Painleve analysis enables us to integrate equations of motion 
for the Goryachev-Chaplygin top W\ as well. It is another integrable system on e*(3) 



at C2 = 0. In this case we have to introduce some seven-dimensional system with the five 
constants of motion. Then this seven- dimensional system may be reduced to the Toda 
system The similar relations between the Toda flow and the integrable system on e*(3) 
is discussed in [jl4|. Now we want to compare the Toda flow with another integrable system 
on e*(3). 

Let us turn to the deformation of the completely symmetric top ( |1.10|) . In 0, the 
following transformation of the independent time variable was proposed 

t^u: ± = ^{y,y)l (3.23) 
au 6 at 

because of the "weak" Kowalewski-Painleve criterion. In the new time variable the initial 
Euler-Poisson equations ( p..l4| ) are 

d ^ 1 d , , 



Namely these equations were integrated in hyperelliptic quadratures in [|T^ . 

If we want to compare integrable system on the sphere S"^ with the Toda system, note 
that equations ( |3.19| ) are invariant by 

1 

Zj I > — , I — ^ —Zj+3 , J — 1; 2, 3 . (3.25) 



8 



The first equation in ( |3.24]) has the similar property by 



Vj ' ) Jj ' ^ Jj ^ J — 1? 2, 3 . 

Using this observations we introduce analog of the Flaschka variables 

Sl = Vi"^ , S2 = , S3 = , Si S2 S3 = 1 

se = yiJi, "54 = ?/3 J3 , S5 = y2 J2 ■ 
satisfying the following equations 



d 

du 


= Sl 


>5 


- -54) , 


d 

du 


2 




-S3) 


+ 


S4 

2 


{S5 


- Sq) 


d 

du 


= S2 


>6 


- S5) , 


d 

du ^ 


2 






+ 


S5 
2 


{S6 


- S5) 


d 

du 


= S3 


>4 


- Se) , 


d 

du 


a' 
2 


is3 


-S2) 


+ 


se 
2 


(S4 


- S5) 



(3.26) 



The first column of equations is completely coincided with the corresponding Toda equations 
( p.l9| ). The second columns differ on the at most quadratic polynomials. Thus, for these 
equations we can directly apply the linearisation procedure due by Adler and van Moerbeke 
^ . The freedom in definition of the variables Sj may be used to make the Laurent solutions 
a bit easier. 

The constants of motion for the flow ( p.26| ) are given by 



S*! — Sl S2 S3 — 1 , 

5*2 = S4 + S5 + S6 = C2 = , 

'S's = ( S4 S2 + S5 S3 + Sg Sl ) - (si S2 + Sl S3 + S2 S3) = ai , (3.27) 
5*4 = S4 S5 Sg + a (si Sg + S2 S4 + S3 S5) = bi , 

The three constants 5*1, 5*2 and 5*4 coincides with the corresponding Toda constants. The 
Hamiltonian 5*3 is a qubic polynomials now. 

There is the permutation invariance similar to the Toda lattice ( |3.21| ). Thus, the invariant 
manifold defined by ( p. 271) possesses the third order automorphism ( |3.22|) , which simplifies 



the linearisation procedure. 

It may of course be quite difficult to find variables similar to {zj} or {sj} associated with 
a given mechanical system. These variables have to satisfy some special conditions 0. 
For instance, the corresponding equations of motion have to include at most second order 
polynomials only,see ( p.l9| , |X^ ). However, if we can introduce such variables, the Adler and 
van Moerbeke methods JI], ^, |16| enable us to integrate a given mechanical system. 

Thus, for a motion on sphere ( |1.13|) one has to embed the afiine invariant surface defined 



by ( p.27| ) into the projective space, whose closure is a principally polarised Abelian surface. 
It enables one to define the system to linearising variables. Then we have to prove that the 
vector field corresponding to 5*4 ( |3.27] ) gives the highest flow with respect to the same hyper- 
elliptic curve of genus two. It will complete the linearisation of the integrable deformation of 
the spherical top. Of course, this general machinery leads to the particular results obtained 
in m. 
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4 Conclusion 



In the algebro-geometric approach due to Adler and van Moerbeke [||, 16|, the algebraic 
curve may be constructed without any Lax pair representation. For the considered motion 
on sphere (|1.1CI| ), by substituting the Laurent solutions in the invariants ( p. 27] ) one gets the 



hyperelliptic curve [|T0|. Starting with this curve and the linearising variables |T^ the 2x2 



Lax pair may be obtained (see [T^ for a review). 

In this note we tried to construct the Lax pair in framework of the group-theoretical 
approach to integrable system ||13| , ^ . Applying the method of finite-band integration to the 



adequate Lax matrix, we hope to get solutions, which may be more simpler than the original 
formulae |T^, as for the Kowalewski top 0]. 



Further properties of the integrable deformation of the completely symmetric top, like 
action-angles variables, Poisson structures of the seven- dimensional system and separation 
of variables are under study. The results and the more detailed geometric description will 
be pubhshed elsewhere. 
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